Form Factors for Semileptonic Decays of Heavy-Light Hadrons 
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The strong coupling lattice QCD solution jl| for the Isgur-Wise functions, pctrametrising the semileptonic decay 
form factors of hadrons containing a single heavy quark, is reviewed. Several useful features connected with the 
result are pointed out. 



In the quark mass Af — > oo limit, QCD has 
an exact SU(2Nf) spin-flavour symmetry. This 
symmetry relates many matrix elements and form 
factors of hadrons containing heavy quarks ^,|| . 
In particular, the semileptonic weak decay form 
factors of s— wave hadrons containing a single in- 
finitely heavy quark can be reduced to two un- 
known functions, one for mesons (^) and another 
for baryons (C). The knowledge of these Isgur- 
Wise functions is of practical importance in a pre- 
cise determination of the quark mixing matrix el- 
ement Vcb- the experimentally measured B D 
semileptonic decay rates can be fitted to a leading 
term involving ^ plus phenomenological estimates 
of corrections suppressed by powers of l/M . 

In the Af — > cxD limit, it is convenient to scale 
all variables so as to explicitly separate the factors 
of heavy quark mass, e.g. a heavy hadron state 
is characterised by its four-velocity and the 
Isgur-Wise functions depend only on the Lorentz 
invariant v ■ v' . In the static geometry, the light 
QCD degrees of freedom decouple from the dy- 
namics of the heavy quark, and flavour indepen- 
dence of QCD fixes the absolute normalisation of 
the form factors: ^(w • w' = 1) = 1 = ({v ■ v' = 1). 

As an explicit case, consider the form factor 
for semileptonic weak decay B — > DlV, and its 
behaviour as m?, and TOc go to infinity. Only the 
partially conserved vector current contributes: 

{D{p')\cj,b\Bip)} = ip + p')J+{q') 



Here q = p — p' is the momentum transfer, and 
vv' = (M| + Ml, - q^)/2MBMD. To eliminate 
the spurious pole at = 0, we must have /o(0) = 
/+(0). The divergence of Eq.(l) gives 



mb - rUc 



M^oo 



(1) 



{D{p')\cb\B{p)) = {MB + MD)fo{q^) 

VMbMi^ {l + v v') £,{v ■ v') . (2) 

This result contains several noteworthy points: 
(a) A partially conserved current can undergo an 
ultraviolet finite renormalisation. Thus in general 
a renormalisation constant Zo{Mb, MB,q^) ap- 
pears multiplying the matrix element on the l.h.s. 
of these equations. According to the Ademollo- 
Gatto theorem the deviation from 1 of the 
form factor for the charge operator is of second 
order in the symmetry breaking parameter: 

f+{MB^ MD,q^ ^0) = l + 0((mb-m,)2) .(3) 

Since only first order deviation from the heavy 
quark flavour symmetric limit appears in Eq.(2), 
it follows that the associated renormalisation con- 
stant obeys Zo{Mb = MD,q^ = 0) = 1. In 
a mass independent renormalisation scheme, e.g. 
dimensional regularisation, the renormalisation 
constants are independent of the momentum in- 
sertion. In such a case, we can forget about Zq 
altogether provided we use Eq.(2) in the heavy 
quark flavour symmetric limit only. On the lat- 
tice, with an explicit momentum cutoff, a q^ 
dependent renormalisation constant is not ruled 
out. But the fact that the quark mass opera- 
tor is a local point operator on the lattice (note 
that Fourier transform of a (5— function is a con- 



stant) may still allow one to "infer" the contin- 
uum result, particularly when the lattice result 
has a simple structure. 

(b) Due to an enhanced kinematic factor on the 
r.h.s., the quark-hadron duality analysis of Eq.(2) 
produces a sum rule yielding a tighter upper 
bound for ^ than Bjorken's ||]: 

iiv-v) < 2/{l + vv') . (4) 

An implication is that the sum rules containing 
the scalar operator would have less contamination 
from excited states and saturate faster compared 
to sum rules containing other operators. 

(c) The 3— point function appearing on the l.h.s. 
obeys a Ward identity, which relates it to 2— point 
functions 

(mfc - m,){D{y)\cb{qf, - 0)\B{x)) 

= {D{y)\D{x)) - {B{y)\B{x)) . (5) 

This identity guarantees that {Mb — Me, , — 
0) = 1. Moreover, it is exact even for distinct nib 
and mc, and therefore contains enough informa- 
tion to determine ^ (even though = 0, w-t;' can 
be varied by changing Mb and Md)- In practice, 
the difference of 2— point functions may be a sim- 
pler object to study than the 3— point function, 
for example while applying QCD sum rules. 

In Ref.||l|, ^ was evaluated in strong couphng 
lattice QCD, and then the result was converted 
to continuum language. The logic of this analysis 
was inspired by Wilson's renormalisation group 
based solution of the Kondo problem — the 
weak and strong coupling results for the same 
quantity can be related to each other provided 
one can evaluate the RG connection between the 
two with high precision. It must be kept in mind 
that the weak and the strong coupling fixed points 
describe quite different physics; the scaling rela- 
tions which hold at one fixed point may not hold 
at the other fixed point. Generically, we can write 

f{\) = r[l + 0{a/X)] , r = /(A-oo) , (6) 

where A is the correlation length. Scaling is ex- 
act only on the critical surface defined by A = cxd. 
The 0{a/\) terms are non- universal, i.e. they 
can be different for different quantities and hence 
violate scaling. As a consequence, even along the 




Figure 1. Strong coupling diagrams for the scalar 
form factor in S ^ DlV decays: (a) the leading 
tree level contribution, (b) an example of the next 
order loop contribution. 

renormalised trajectory (i.e. for an action con- 
taining no lattice artifacts, yet having a finite 
cutoff), there exist scaling violations. For exam- 
ple, Wilson's calculation demonstrates how the 
spectrum of the Kondo problem changes under 
RG evolution. (Another example is the ratio of 
two particle binding energy to a single particle 
mass. It vanishes at the strong coupling fixed 
point, while it can be non-zero at the weak cou- 
pling fixed point of an interacting theory.) Thus 
scaling relations of the theory must be extracted 
only in the weak coupling region, after individu- 
ally converting each quantity calculated at strong 
coupling to its weak coupling analogue. 

In QCD scaling violations are ©(oAqcu), once 
the coupling is traded off for a dynamically 
generated scale. To keep them under control, it 
is of paramount importance to select a quantity 
which has a simple RG connection between the 
strong and weak coupling fixed points, and which 
still contains the physics of interest. The choice of 
the quark mass operator happens to be crucial for 
this reason — in a lattice regularisation respecting 
the quark flavour symmetry, it is a local point 
operator with a non-perturbatively constrained 
renormalisation constant. 

The scalar form factor is easily evaluated in 
the strong coupling limit, the leading term be- 
ing a simple pole contribution from the tree level 
hadron diagram of Fig. la. The sublcading correc- 



tions arise from diagrams containing light quark 
loops. The diagrams where the loops are part of 
the external B and D legs merely redefine the ex- 
ternal parameters. The modification of the form 
factor arises from diagrams of type Fig. lb, where 
the loops interact with the scalar hadron state. 
All such diagrams involve at least two more heavy 
quark propagators compared to the tree level di- 
agram, and hence are suppressed in the M oo 
limit given the mass independence of quark-gluon 
interactions. 

Explicitly, with staggered fermions (their chiral 
symmetries are essential in restricting the renor- 
malisation of the quark mass operator) in the 
strong coupling limit, the scalar form factor is: 



1 



(7) 



where Msc is the mass of the scalar meson cb. 
The normalisation of this form factor at = 
confirms the expected non-renormalisation con- 
straint. The "sin" —function in the denomina- 
tor is a reflection of the nearest neighbour lattice 
action used in the analysis; a more general lat- 
tice action would give rise to a different function. 
(As a matter of fact the explicit structure of the 
strong coupling lattice action generated by RG 
evolution is not known for QCD.) What would re- 
main unchanged for any lattice action, however, 
is the simple pole structure of fo{q^)- In the con- 
tinuum language, this simple pole corresponds to 



(8) 



where we have set Mb ~ Mr, and then taken the 
M ^ oo hmit (which implies Msc —> Mb + Md)- 
A similar analysis for the baryon form factor 
yields C,[v ■ v') = 2/{l + v v'). 

The result of Eq. (8) fits the experimental data 
for B D*lV decays reasonably well ||,^. It 
is reassuring to find that the agreement is within 
~ 10%, which is the precision of the experimen- 
tal data as well as the estimated magnitude (e.g. 
from QCD sum rule calculations) of the 0{1/M) 
difference between the actual form factor and ^. 

Encouraged by this agreement, we turn the 
problem around and construct a phenomcnologi- 



cal model of hadron wavefunctions which repro- 
duces Eq.(8) exactly. Such wavefunctions can 
be useful for studying other hadronic properties. 
The scalar form factor at = is nothing but 
the overlap of the initial and final state light-cone 
wavefunctions. Eqs.(2) and (8) require 

■ / 2 
du 4>Uu) ct^Biu) = — ; . , (9) 



where is the fraction of longitudinal light- 
cone momentum carried by the heavy quark. A 
simple choice for the wavefunction is 

0(m) = 2{M/mof'^ u e""*^/™" , (10) 

with an Af— independent parameter rriQ. By re- 
placing the factor of u with 1 — e~" and choosing 
mo to be two-thirds the constituent quark mass 
in the chiral limit, the structure of this wavefunc- 
tion can be symmetrised between the two con- 
stituent quarks of the meson. It would be inter- 
esting to explore how accurately this wavefunc- 
tion describes many other properties of hadrons. 
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